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Bistability generated via a pure noise-induced phase transition is reexamined from the view of 
bifurcations in macroscopic cumulant dynamics. It allows an analytical study of the phase diagram in 
more general cases than previous methods. In addition using this approach we investigate spatially- 
extended systems with two degrees of freedom per site. For this system, the analytic solution of the 
stationary Fokker-Planck equation is not available and a standard mean field approach cannot be 
used to find noise induced phase transitions. A new approach based on cumulant dynamics predicts 
a noise-induced phase transition through a Hopf bifurcation leading to a macroscopic limit cycle 
motion, which is confirmed by numerical simulation. 
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INTRODUCTION 



The interplay between nonlinear dynamics and noise 
often generates interesting and counterintuitive phenom- 
ena. Popular examples are stochastic resonance co- 
herence resonance |2j, and noise induced phase transi- 
tions 0, 0, . In the latter example noise creates an or- 
dered phase which does not exist in the absence of noise. 
Unlike noise induced transitions in systems with few 
degrees of freedom, the noise induced phase transition 
breaks ergodicity and has the characteristics of a genuine 
phase transition. In previous studies, many variations 
of pure noise-induced phase transitions were introduced 
0. Spatial patterns can be induced via the pure noise- 
induced phase transition |E S S • The noise- induced 
first or der p hase transition was also shown to be possible 
[lOl IllL Il2j . The systems with colored noise were inves- 
tigated by various groups Furthermore, 
the bistability created by the noise-induced phase transi- 
tion exhibits stochastic resonance when a time-periodic 
perturbation is added UM a nd it can lead to propaga- 
tion of harmonic signals [l7| . The idea of noise-induced 
phase transition was also used in coupled Brownain mo- 
tors, mini 



Most of previous investigations take a mean field ap- 
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proach and use the self-consistent condition 0, l2pj | 

/CO 
xP st (x;(x))dx. (1) 
-CO 

to determine the mean (x) and also bifurcation points. 
This method yields an exact solution for the phase 
boundaries. However, to solve Eq. one needs to know 
the analytic solution of the stationary Fokker-Planck 
equation P st (x;{x)), which is in general not available. 
Furthermore, if the system does not have a stationary 
state (i. e., (x) is time-dependent), we must use time- 
dependent self-consistent condition which is prohibitively 
more difficult. For certain types of spatially extended 
stochastic problems, there is an approximate method 
which replaces stochastic dynamics with effective deter- 
ministic dynamics [2l[. However, the extent of appli- 
cability to other models is not known. A general and 
systematic method is highly desired. 

In this paper we present a systematic method to inves- 
tigate noise-induced phase transitions. While it does not 
provide exact solutions, the method does not require an 
analytical expression of the stationary state probability 
distribution and can be applied to general cases includ- 
ing time-dependent problems. In the following section, 
we investigate a model system with a single variable in- 
troduced by Van den Broeck et al. 3] for which an exact 
solution is known. The present method predicts a pitch- 
fork bifurcation to an ordered state as the noise intensity 
increases and also the reentrant transition to a disordered 
phase at a higher noise intensity. This behavior is quali- 
tatively in a good agreement with the exact results. 

Then, we apply the same method to a model with two 



2 



variables. The model is expected to undergo a noise- 
induced phase transition to a time-dependent ordered 
phase for which the time-dependent self-consistent ap- 
proach is not practical. The present method predicts a 
Hopf bifurcation to a macroscopic limit cycle phase from 
a disordered state as the noise intensity increases and 
shows also a reentrant transition. We also demonstrate 
that the present method can provide other information 
such as a period and amplitude of the oscillation. 



II. NOISE-INDUCED PITCHFORK 
BIFURCATION 

In this section, we consider the following stochastic 
system of N globally coupled microscopic variables {xi}: 



*i = f(xi) ~ j- ~ x i) + 9^i)iiif) (2) 

where D is a coupling strength and is a Gaussian 
white noise, defined by 

&(*)> = o, UtMt')) = * 2 SiAt - 1')- (3) 

Equation (0) is interpreted in the Stratonovich sense. 

With certain nonlinear functions f(x) and g(x), Eq. (0 
exhibits a phase transition from a disordered phase 
((x) = 0) to an ordered phase ((x) ^ 0) as the noise inten- 
sity a 2 increases. At larger noise intensities the system 
undergoes a reentrant transition to another disordered 
phase ((x) =0) |3,E||- Since the ordered phase does not 
exist in the absence of noise, the phenomena is called pure 
noise-induced phase transition. Van den Broeck et al. 
|3| originally investigated it using the following nonlinear 
functions: 

f(x) = -x(l + x 2 ) 2 , g(x) = l+x 2 (4) 

because a stationary state solution to the correspond- 
ing Fokker-Planck equation can be obtained analytically, 
which allowed them to find the exact phase boundary 
with the self-consistent equation QJ. 

In the following, we present a method which allows us 
to investigate more general cases approximately but with- 
out an analytical probability distribution. We apply the 
method to this model Q and compare the results with 
the exact solution and also with numerical simulations. 



A. Moment Dynamics 

Assuming TV — > oo, we write Eq. in a mean field 
form: 

x = f(x)-D(x-(x))+g(x)at)- (5) 



Taking the mean of Eq. © under the Stratonovich inter- 
pretation, the dynamics of (x) is given by 

(x) = (f(x)) + ^(g'(x)g(x)). (6) 

Expanding f(x) and g(x) in Taylor's series around (x), 
Eq. © forms an infinite set of simultaneous ordinary 
differential equations: 

& = E S {/ (n) «*» + y [s'((*>M(*>)] (n) ) (7) 

71 = ' ^ ' 

(l n -Dnfln 

+ E {/ (m) ^» + yb'(<z>M<*>)] (m) } 

+ n( "" 1 ^ yb 2 (W)l' m) . (8) 

Here, is the n-th order derivative and ji n = ((x — 
(x)) n ) the n-th central moment with, by definition, /iq = 
1 and ji\ = 0. 

Some previous studies 0, have considered only the 
0-th order term, thereby neglecting the fluctuations x — 
(x), and worked with the equation 

(i) = /«*» + yj'(WMW). (9) 

Since this equation does not depend on the coupling con- 
stant D, it cannot explain the reentrant transition. In 
fact, Eq. (0 is exact when D — > oo under which condi- 
tion no reentrant transition takes place. |3|,l| For a finite 
D, one must consider higher order terms. 

In order to see how the higher order terms create the 
reentrant transition, we investigate the previous model 
Eq. (0J). Here we show the equations of motion only for 
(x) and /i 2 : 

(x) = (a 2 - 2)/i 3 - /i 5 

+ [a 2 - 1 + 3(a 2 - 2)/i 2 - 5/i 4 ] (x) 

- 10/i 3 (x) 2 + (^ 2 -2-10/i 2 )(x) 3 - (x) 5 (10) 

/i 2 = (J 2 + 2[2cr 2 - (1 + D)]fi 2 + (3a 2 - 4) /i 4 - 2/i 6 
+ [2(5a 2 - 6)/i 3 - 10/i 5 ] (x) 

+ [2(J 2 + 12(a 2 - l)/i 2 - 20/X 4 ] (xf 

- 20/i 3 (x) 3 + (a 2 -10/i 2 )(x) 4 . (11) 

Even with higher order terms, Eq. (^UJ) still does not 
explicitly depend on the coupling constant. The effect 
of coupling arises from the dynamics of the second and 
higher order moments as Eq. (|11|) indicates. 

Noting f(—x) = —f(x) and g(—x) = g(x), the system 
(O is invariant under the variable transformation x — > 
—x. Due to this symmetry the odd moments must be zero 



FIG. 1: Stationary moments (left panel) and \i\ (right 
panel) obtained by Gaussian approximation (solid lines) and 
numerical simulation (dashed lines). The coupling strength 
D = 10 is used. The Gaussian approximation appears in a 
good agreement with the simulation below a 2 = 2. However, 
it overestimates both 2nd and 4th order moments above a 2 = 
2. This rapid growth causes reentrant transition at a smaller 
noise intensity than the exact solution. The grey region shows 
the location of the ordered phase. 
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FIG. 2: The phase diagram. The ordered phase ((x) ^ 0) 
is above the lines and the disordered phase ((x) = 0) below 
them. The dashed line is the exact solution obtained by Eq. 
Q] the solid line is the Gaussian approximation. 



when (x) = 0. Then, we find a fixed point at (x)* = 
and /i2n+i = 0- The even moments at this fixed point 
are not zero and must be determined by setting the r.h.s. 
of Eq. ((HJ) to zero. For example, Eq. (fTTJl provides the 
following equation: 

a 2 + 2[2a 2 - (1 + 2?)]/x* + (3a 2 - 4) /4 - 2/i*. = (12) 

A linear stability analysis of Eq. ((!())) yields the bifur- 
cation condition for a pitchfork bifurcation: 

a 2 c - 1 + 3(a 2 c - 2)£ -5/4 = (13) 

where a c is a critical noise intensity. Even without the 
exact knowledge of the higher moments useful informa- 
tion can be derived from Eq. (|T3|). Since both \±2 and {14 
are non- negative, l.h.s. is alway negative for a 2 < 1. In 
this regime the fixed point (x)* =0 is stable regardless 
of the magnitude of D. For 1 < a 2 < 2, only the 0-th 
moment term is positive but the other terms are nega- 
tive and support stability. The bifurcation is possible in 
this range of noise intensity only when the moments are 
sufficiently small. Since increasing the coupling strength 
reduces fluctuation, the bifurcation takes place above a 
certain magnitude of D. 

Interestingly, the role of the second moment term 
changes at a 2 = 2. At higher values it supports instabil- 
ity of the fixed point. The 4-th moment term is always 
negative and supports stability. When n± grows faster 
than 112 with increasing a 2 , it eventually dominates and 
l.h.s. of Eq. ((13)) becomes negative again. Then, the sys- 
tem reenters the disordered phase. In order to determine 
the critical noise intensity a c from Eq. ((13)). one needs 
to know stationary even moments \i\ and which are 
to be determined by Eq. ((12)) . In turn, it requires \±*§. 
At the end, all even moments must be simultaneously 



solved, which is practically impossible for general cases. 
An approximation is necessary. 



B. Gaussian Approximation 

When f(x) or g(x) is nonlinear, the stochastic dynam- 
ics (0 is not a Gaussian process and in general we can- 
not solve the system of equations and JSJ) exactly. 
For an approximate solution it is convenient to assume 
a probability distribution P(x; (x)) for which the cumu- 
lants above a certain order are negligible. We chose here 
the simplest example where the cumulants above the sec- 
ond order are set to zero (Gaussian approximation). In 
general, this approximation is quantitatively not justified 
but reproduces the main features of the noise-induced 
phase transition, especially the reentrant transition into 
the disordered phase. 

In the Gaussian approximation all odd moments van- 
ish. The even moments can be expressed by the 2nd 
moment fi2- For example, 

/i4 = 3/i?, an d /i6 = 15/^2- (14) 
Applying these relations, Eq. ([T2|) becomes 

a 2 + 2[2(j 2 - (1 + D)]H2 + 3 (3a 2 - 4) (/i*) 2 - 30(/i*) 3 = 

(15) 

which determines function of a 2 . In turn, we can 

determine higher order even moments via the Gaussian 
approximation. Figure Q plots \i\ and /i| obtained by 
the Gaussian approximation and also the results of nu- 
merical simulation for comparison. Up to a 2 = 2, the 
Gaussian approximation is in a good agreement with the 
simulation. At higher noise intensities, non- Gaussian be- 
havior becomes large and the Gaussian approximation 
significantly overestimates the fluctuation. 
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FIG. 3: The mean (x) by the Gaussian approximation (solid 
line) and by numerical simulation (dashed line). The coupling 
strength D=10 is used. 




FIG. 4: Phase diagram for the two variable case. In the region 
above the lines we observe oscillations of the mean, below it 
we do not. The line labeled theory is the analytic solution 
in the Gaussian approximation, the line labeled simulation is 
the numerically obtained result for the Langevin dynamics of 

Eq.[T7|and[THlfor N 625 s y stems - 



Once we find all stationary moments, we can quantita- 
tively evaluate the stability condition (p"3|) for (x). With 
the Gaussian approximation (|14[) . the bifurcation condi- 
tion (|T3|) becomes: 

a 2 c - 1 + 3(<x c 2 - 2) M * - I5(p* 2 ) 2 = (16) 

Here, \i\ must also satisfy Eq. ([T5|) with noise intensity 
<j c . In other words, Eqs. (fT5|) and (|16|) must be solved 
simultaneously for a c and /ij. 

Figure shows the results and compares them with the 
exact solution obtained from Eq. Although there is 
a clear quantitative discrepancy between the approxima- 
tion and the exact solution the main features, namely the 
entrant transition into the ordered phase at medium noise 
intensities and the reentrant into the disordered phase at 
high intensities is reproduced. As in the exact solution 
a certain minimum coupling strength is required for the 
transition to take place. 

We can also evaluate the mean (x) as a function of a 2 
from Eqs. 0, OH), and 0. The result is plotted in 
Fig. [3| along with the results of simulation. The agree- 
ment near the first transition point is quite good. How- 
ever, the Gaussian approximation predicts the reentrant 
transition much earlier than the simulation as mentioned 
earlier. 

We noticed that the present results resemble the phase 
boundary for the 2-dimensional system with a local cou- 
pling examined in Ref. surprisingly well. This coin- 
cidence is partly due to the fact that the locally coupled 
system has much larger fluctuation than the globally cou- 
pled system, which induces the reentrant transition at a 
smaller noise intensity. Since the Gaussian approxima- 
tion overestimates the fluctuation, it shares some simi- 
larity with locally coupled systems. 



III. NOISE-INDUCED LIMIT CYCLE 

In this section we investigate a model with two vari- 
ables at each site. Multivariate stochastic systems are 
mathematically quite difficult. Even if the system has 
a stationary state, it is usually hard to find an analyti- 
cal expression of its probability distribution. It will be 
even more difficult if the system does not have a station- 
ary state and the probability distribution is explicitly 
time-dependent. The lack of an analytical expression of 
the probability distribution makes the standard method 
based on the self-consistent equation (Q futile. Direct 
numerical simulation of multivariate Langevin equations 
demands more computational power than the single vari- 
able cases. However, we expect that the Gaussian ap- 
proximation provides a degree of accuracy similar to that 
in the single variable case without increasing mathemat- 
ical difficulty. 

We use the following simple model that keeps a close 
connection to the previous model: 

Xi = f(xi) - ^ ^2(xi - xj) + g(xi)£i(t) - yi (17) 

3 

Vi = kxi (18) 

with f(x) and g(x) defined by Eq. (@J) as before. Again, 
we expand the dynamical equations in terms of the cen- 
tral moments: 

<*> = E^r{/ ( " ) (( a; )) + Y[5 , (( a; ))5((^))] ( " ) } 

n=0 ' ^ ' 

- (y) (19) 



(y) = Hx) (20) 
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£2,o = <x 2 + 2 [2<x 2 - (1 + £>)] /i 2 ,o 

- 2/Ji.i + (3a 2 - 4)/t 4)0 

- 2/t 6 ,o + 2 [(5<j 2 - 6)/t 3) o - 5/t 5) o] (a;) 
+ 2 [o- 2 + 6(<r 2 - l)/t 2 ,o - 10au,o] (^) 2 

- 20/i 3) o(a;} 3 + ((T 2 -10/i2)(a ; } 4 



(23) 






£0,2 = 2fc/ti,i 

Mi.i = (<t 2 - -D - l)/fi,i + fc M2,o - Mo,2 + (<? 2 

- M5,l + [3(cr 2 - 2)/i 2 ,i - 5/i 4 ,i] (x) 
+ [3(a 2 - 2)a*i,i - 10M3,i] {*? 

— 10/t2,i (a;) 3 - 5/ti,i(a;) 4 . 



(24) 



2)^3,1 



(25) 



Here only the lowest order moments are shown. 

Taking into account the symmetry of the model, there 
is at least one fixed point at (x)* = (y)* = and 
Mn 2m+i = M2m+i n = 0. Stationary even moments are 
determined by an infinite set of simultaneous equation. 
Here we show three conditions derived from Eqs. ([25])- 
C3: 



a 2 + 2 [2a 2 - (1 + £>)] /i* + (3a 2 



4)m1,o - 2 Me,o = 0. 

(26) 



FIG. 5: Snapshots of an ensemble of 625 elements in phase- 
space. The single dots denote the individual elements. The 
solid circle shows the mean and the line depicts the trajectory 
of the mean. The time series goes from left to right and from 
the top to the bottom. Parameter values D = 10 and k = 0.1 
are used. 



An,m — TiD fJj n ,m ^n-l,m+l 

+ rafc/i n+ i >m _i + mkfi ni m-i(x) - n/i n _i >m (2/) 

£ {/«>«*» + 

(21) 



^(n-l)to- 2 , m ^ 2((i))|(t) 



where n+ra > 2 and /i n , m = ((x— < x >) n (y- < y >) m ) 
with /io,o = 1 an d Mi,o = Mo,i = 0. 

For the present model Eqs. ([T5)l and ([2T)l are ex- 
plicitly written as 



= (°" 2 - 2)/*3,o - M5,0 

+ [cr 2 - 1 + 3(<T 2 - 2)/l 2 ,0 ~ 5/i 4 ,o] (x) 

- 10/i 3 ,o(a;) 2 + (a 2 -2-10/i2,o)(a;) 3 

- (x) 5 -(y) (22) 



2,0 



MO, 2 



+ (a 2 



^3,1-^5,1 



(27) 



0. 



and n\ 

Now, we apply the Gaussian approximation to this sys- 
tem. For any two variable Gaussian system, all odd mo- 
ments are zero (/i m , n = for m + n = odd integer) and 
any even oder moment can be expressed as a product 
of the 2nd order moments, /i2,o> Mo, 2, and fii^. For the 
present model we need only the following relations: 



M3,l — 3^2,0/^1,1? M4,0 - 
M5,l = 15/i2 5 oMl,l 5 M6,0 



3/^2,0 

= 15/i| . 



(28) 



Under the Gaussian approximation, the stationary even 
moments o an d Mo 2 can be determined by Eqs. ([25)1 
and ([27)1 . 

The equations of motion (|22|) - (|25|) be- 
come a five-dimensional dynamical system of 
(?/) /i 2 ,o Mo, 2 Mi,i)- A standard linear stability 
analysis yields a Jacobian: 



/Ai -1 \ 

fc * 

A 2 

2k 

\0 fc -1 A 3 / 



(29) 



where 

Ai = 
A 2 = 
A 3 = 



2(2<r 2 



- 3(<x 2 - 2) M ^ - 15(/4 ) 2 (30) 
1 - £>) + 6(2<x 2 - 4)^ >0 - 90(/4, $l) 
-D + 3(cr 2 -2)^ -15(/i*2,0) 2 . (32) 



6 




FIG. 6: Power spectrum of a limit cycle motion obtained by 
numerical simulation (solid line) and the period obtained from 
the linear stability analysis (dashed line). Parameter values 
are D=10, cr 2 =2, and k=0.1. The agreement is perfect. 



The Jacobian ([25)1 is in a block diagonal form and the sta- 
bility of (x) and (y) are separated from that of the higher 
order moments, which makes analytical investigation eas- 
ier. The 2-by-2 block at the top- left corner determines 
the stability of (x)* = (y)* =0 and its eigenvalues are 
given by 

A = 1 (Ai ± \J Af - 4k) , (33) 

which indicates that the fixed point becomes unstable 
at Ai = 0. This bifurcation condition is identical to 
Eq. ([TB)l and therefore, the two- variable model also un- 
dergoes reentrant transition because of the same reason 
as in the single variable model. Figure 0] compares the 
phase boundary obtained by the Gaussian approximation 
and the simulation results. Quantitatively, the disagree- 
ment is rather large. However, the qualitative features 
are correctly captured. 

Since the eigenvalue has an imaginary part at the bi- 
furcation point, it is a Hopf bifurcation and a stable Limit 
cycle is formed above the critical noise intensity, which 
is confirmed by numerical simulation. Figure shows 
snapshots of an ensemble of particles. One observes a 
quite regular limit cycle motion of the mean, even though 
the system is rather small (N = 625) and the individual 

I 



units are spread widely around the mean. The width of 
the cloud in x direction is much larger than the one in 
y direction due to small k, which is in a good agreement 
with Eq. J2ZJ). Analogous to the single- variable case we 
will call the oscillating phase the ordered phase and the 
non-oscillating phase the disordered phase. 

From the imaginary part of the eigenvalue the period 
of oscillation near the bifurcation point is approximately 
given by T = 2tt j\[k. Since it does not depend on any 
moment, this period is valid even without the Gaussian 
approximation. Indeed it perfectly agrees with numerical 
simulation as shown in Fig. El 

If the time evolution of the means is needed, Langevin 
equation or time-dependent Fokker-Planck equation are 
usually solved numerically. Numerical simulation of cou- 
pled Langevin equations is computationally rather time 
consuming especially near the bifurcation point due to 
the finite size effect. One needs a large number of sam- 
plings to obtain a reasonable statistics. Furthermore, 
ensemble averaging is cumbersome since each realization 
oscillates in a different phase. Numerical integration of 
time-dependent Fokker-Planck equations does not have a 
problem of statistical error but often suffers from numer- 
ical instability. A special care may be needed. While it is 
difficult to solve moment dynamics (|22|). ([20)) and ([23)) - 
([25)) analytically even with Gaussian approximation, it is 
much easier to solve them numerically compared to the 
Langevin or Fokker-Planck equation. Since the moment 
dynamics is deterministic, there is no statistical error. 
The upper panel of Fig.[7|shows the time-evolution of (x), 
(y), and /i2,o with the Gaussian approximation. Other 
moments, /io,2 and /ii 5 i (not shown) converge to station- 
ary values. The lower panel shows the results of numeri- 
cal simulation. Only one realization with 10000 particles 
is shown. Although the amplitude and period are over- 
estimated by the Gaussian approximation, all qualitative 
features are captured. 

Finally, we discuss a special case where (relax- 
ation oscillation limit). In this limit, Eq. (|20j) indicates 
(y) varies very slowly. Furthermore, from Eq. (|27|) the 
fluctuation of the variable y is negligibly small, suggest- 
ing that all particles experience the same value of y. Since 
x dynamics is much faster than a "stationary" prob- 
ability distribution is formed before y varies. In another 
word, the variable y in Eq. JTJJ) is just a parameter for 
the dynamics of x. In this case, we can investigate the 
Hopf bifurcation using the self-consistent equation QJ. 
Following Ref. |3| , the stationary distribution is given by 



p ^ {y)) = z exp L dx \ x ~ TTT^ - (i + oV J 
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FIG. 7: Time evolution of (x) (solid line) and (y) (dashed 
line) obtained by numerically solving (|2()f) and (|22^ - (|25| in the 
Gaussian approximation (upper panel) and Langevin equa- 
tions |T7)) and (fT8|) (lower panel). For the Gaussian approx- 
imation /i2,o (thin line) is also shown. Parameters D = 10, 
a 2 = 2 and = 0.1 are used. 



where Z is a normalization constant. The self-consistent equation (QJ yields a relation between (x) and (y): 



1 / i+0 ° 
(x) = — J dx- 



1 + ar 



■ exp 



L> 1 

(T 2 <J 2 1 + X 2 



D(. 



x) - jy) f x 
a 2 \ 1 + a 



arctan x 



(35) 



which corresponds to a nullcline of the relaxation oscil- 
lation. Figure [HI plots the nullcline and an actual tra- 
jectory from numerical simulation. The trajectory fol- 
lows one branch of the nullcline and jumps to the other 
branch analogous to the deterministic relaxation oscil- 
lation. This nullcline is exact. However, its derivation 
requires an analytical expression of the probability dis- 
tribution and may not be applicable to other cases. An 
alternative method such as the Gaussian approximation 
is useful for general cases where the analytical distribu- 
tion is not available. 



IV. CONCLUSIONS 

We have described a general theoretical method use- 
ful for the investigation of nonlinear stochastic systems 
exhibiting noise-induced phase transitions. It is based 
on the dynamical equations for the central moments. In 
general the exact solution involves the solution of the infi- 
nite dimensional set of ordinary differential equations for 
these moments. However, even without solving this sys- 
tem useful qualitative information can be gathered from 
the equations. Furthermore, quantitatively reasonable 
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solutions can be obtained by neglecting the cumulants of 
the distributions above a certain order and then solving 
the remaining finite set of equations. 

Using this method we investigated two systems. One 
exhibits a noise induced pitchfork bifurcation, the other 
one a Hopf bifurcation. In the latter case, the macro- 
scopic quantities oscillate in time when the system is in 
an ordered phase. This oscillation is purely induced by 
noise via spontaneous symmetry breaking. The macro- 
scopic oscillation suggests a strong synchronization of mi- 
croscopic degrees of freedom despite the presence of noise. 
Actually, it is the noise that generates the macroscopic 
order. On the other hand a strong noise destroys the 
order again. The reentrance into the disordered phase is 
due to the fourth moment, that grows faster with noise 
intensity than the second moment. 

In the Gaussian approximation we have reproduced 
the basic features of the noise-induced phase transition, 
namely the existence of a critical coupling strength and 
the disorder-order-disorder transition. 



Sonderforschungsbereich 555 "Komplexe Nichtlineare 
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FIG. 8: A noise-induced relaxation oscillation. Dashed line: 
a nullcline determined by Eq. (|35f) . Solid line: a limit cy- 
cle trajectory obtained by numerical simulation. Parameter 
values are D—10, a 2 — 2, and k=0.1. The simulation result 
follows the theoretical nullcline and jumps to another branch 
of the nullcline. 
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